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ABSTRACT 

The  prediction  of  torsional  buckling  of  stiffeners  is  a  necessary 
design  skill  in  the  planning  of  any  framed  structure  or  vehicle.   The 
likelihood  of  torsional  failure  is  increased  when  the  stiffener  is 
simultaneously  loaded  from  perpendicular  directions,  a  common  loading 
condition  for  the  frames  of  a  ship  at  sea.   This  thesis  examines  first 
the  case  of  a  single  lateral  load  on  I  beam  and  flat  bar  stiffeners 
and  then  the  instance  of  combined  lateral  and  axial  loadings.   An 
energy  method  is  used  and  formulated  are  expressions  for  the  critical 
buckling  load  for  a  single  laterally  applied  force  and  a  formula 


relating  lateral  and  axial  loads  and  giving  an  indication  of  the 
relative  weakening  of  a  stiffener  caused  by  simultaneous  application 
of  both  forces. 


Thesis  Supervisor:      Dr.  P.C.  Xirouchakis 
Thesis  Reader:  Dr.  J.H.  Williams,  Jr. 
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Length  of  stiffener 

Breadth  of  I  beam  flange 

Angle  of  rotation  of  stiffener  in  the  YZ  plane 

Amplitude  of  cosine  function  for  B 

Height  of  stiffener  shear  center  above  toe 

Warping  Constant 

Additional  lowering  of  the  load  P  due  to  the  load 

being  applied  a  distance  h/2  above  center 

of  the  web 

Young's  Modulus  =  2.068  x  10 8  KN/m 

Shear  Modulus  =  E/2.6 

Height  of  stiffener 

Moment  of  inertia  of  stiffener  about  web  plane 

Local  rectangular  coordinates  at  any  point  x 

on  a  stiffener 

St.  Venant's  torsion  constant 

Dimensionless  parameter 

Moment  about  the  I    axis 

Lateral  load 

Critical  lateral  buckling  load 

Axial  end  load 

Euler  buckling  load 


tf  Flange  thickness 

tw  Web  thickness 

u,v,w  Displacements  in  the  x,y,  and  z  directions 

U  Strain  energy 

W  Work 

x,y,z  Rectangular  coordinates  originating  at 

centroid  of  stiffener 

If  Numerical  factor 

C  Membrane  strain 

v  Poisson's  Ratio 
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INTRODUCTION 

As  is  the  case  with  most  large  engineering  structures,  a  ship's 
lull  and  its  interior  dividing  bulkheads  and  decks  are  basically 
framework  patterns  of  stiff eners  covered  by  an  outer  shell.   It  is 
principally  this  lattice  of  stiffeners  which  determines  the  strength 
and  structural  stability  of  a  vessel.   Structural  soundness  and 
integrity  are  particularly  critical  for  a  ship  because  this 
engineering  structure  must  operate  in  what  can  be  the  harshest  natural 
nvironment  on  earth.   Very  seldom  are  the  frames  of  a  ship's  hull 
ubject  to  loading  in  a  single  plane;  rather,  it  is  most  often  a 
combination  of  loads  caused  by  the  interaction  of  the  various  natural 

nd  man-made  forces  from  different  directions  which  determine  the 
total  stress  on  a  member.   The  subject  of  this  thesis  is  to  develop  an 
expression  relating  the  interaction  of  simultaneous  lateral  and  axial 
Loading  of  a  stiffener. 

Two  of  the  more  commonly  found  shipboard  stiffeners,  the  I  beam 
and  flat  bar  type,  will  be  examined  first  under  the  action  of  a 
centrally  applied  lateral  load  and  then  under  the  combined  actions  of 
:he  central  lateral  load  and  a  per pind icular ly  applied  end  axial  load. 
Concentrated  point  loads  were  chosen  because  they  represent  the  most 
.imiting  and  potentially  dangerous  cases.   An  energy  approach  will  be 
ised  in  all  instances.   First  the  strain  energy  of  the  system  will  be 
ound  and  then  the  work  done  by  the  loading  forces.   By  conservation 
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of  energy,  the  variation  of  the  strain  energy  will  be  equated  to  the 
work  done.   From  the  resulting  expression  will  be  developed,  in  the 
case  of  the  single  lateral  load,  a  formula  for  the  critical  buckling 
load  and,  in  the  case  of  combined  loading,  a  formula  relating  the 
interaction  of  the  lateral  and  axial  forces. 

It  is  hoped  that  study  into  the  problem  of  stiffener  failure 
under  more  realistic  loading  conditions  such  as  this  may  aid  in  the 
formulation  and  upgrading  of  design  codes  used  in  civilian  and 
military  ship  design,  and  any  efforts  that  could  make  a  ship  safer  or 
lighter  would  certainly  be  worthwhile.   To  this  end  it  is  hoped  that 
the  ultimate  contribution  of  this  thesis  is  a  reasonably  accurate 
method  of  predicting,  given  a  value  of  axial  load,  how  much  lateral 
load  may  be  sustained  by  a  stiffener  without  buckling,  or  vice  versa. 
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CHAPTER  1 
DESCRIPTION  OF  MODELS 


The  two  types  of  strengthening  members  that  will  be  considered 
are  the  I  beam  type  and  flat  bar  type,  to  be  known  respectively  as 
Model  I  and  Model  II.   The  coordinate  system  to  be  employed  with  both 
models  is  shown  in  figure  (1-1). 


Z 


o 


->  X 


T77~r7~rrr? 


z/7/7/7 // 


Figure  (1-1)  Coordinate  System 


The  origin  of  the  coordinate  system  will  be  taken  at  the  center 
of  the  beam,  with  u,vf  and  w  being  the  displacements  along  the  x,  y, 
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and  z  axes.   The  deflections  will  be  positive  in  the  direction  of  the 
axes  shown  in  figure  (1-1). 


Torsional  buckling  is  that  failure  mode  characterized  by  a 
deflection  and  twisting  of  the  stiffener  in  the  YZ  plane.   This  type 
of  deformation  is  shown  in  figure  (1-2). 

Y 

As 


»  X 


Figure  (1-2)  Characteristics  of  Torsional  Buckling 
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1.1      Model    I 


Model    I   as    stated    has    been    chosen    to    be    a    simply    supported    I 
beam.      The    case    of    a    beam    with    symmetric    flanges    will    be    considered, 
and    a    lateral    cross    section    of    Model    I    is    shown    in    figure    (1-3). 


fcf 


T" 


w 


±k- 


Figure    (1-3)       Model    I    Lateral    Cross    Section 


The    properties    of    the    Model    I  cross    section    are    listed    below 


14 


c  -  Height  of  shear  center  above  toe  (Flange 

thickness  t,f  considered  negligible  compared 
to  web  height  h) 


c  = 


Cw    -    Warping    constant 


Cw    = 


2      3 
t      \\      b 
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J    -    St.    Venant's    torsion    constant 


J    =  -^-(2btf3  +      htw3) 


Iz  -  Moment  of  inertia  about  the  web  plane 
(contribution  of  top  and  bottom 
flanges  only) 


Iz  =  2  ( ) 

12 
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1.2.  Model  II 


Model  II  is  the  flat  bar  type  stiffener  as  shown  in  figures  (1-4 
and  (1-5). 


Figure  (1-4)   Model  II  -  Flat  Bar  Stiffener 


— » 


<—  tw 


} 


Figure  (1-5)   Model  II  Lateral  Cross  Section 
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The  properties  of  the  flat  bar  cross  section  are  as  follows: 


c  -  Height  of  shear  center  above  toe 


h 
c  =  — 
2 


Cw  -  Warping  Constant 


Cw  =  0 


J  -  St.  Venant's  torsion  constant 


htw 
j  =  


I   -  Moment  of  inertia  about  the  web  plane 


htw 
Iz  =  


12 
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CHAPTER    2 


Model    I     -    TORSIONAL    BUCKLIN3    UNDER    POINT    LATERAL    LOADING 


2.1      Development    of    the    Strain   Energy   Equation 


The    case    of    a    simply    supported    I   beam   of    length    a    subjected    to    a 
point    lateral    load    P    in    the    center    is    illustrated    in    figure    (2-1) 


>  X 


Figure  (2-1)  I  Beam  Under  Point  Lateral  Load 


The  general  equation  for  the  strain  energy  of  a  beam  of  length  a  under 
bending  and  torsion  is  given  by  (refernce  1): 


U  = 


a/2   r— 


■.,4i!."e''»»'p»!  ">» 


G 


Considering  only  that  portion  of  the  beam  to  the  right  of  an  arbitrary 
point  x,  figure  (2-2),  the  forces  acting  will  be  a  single  one  P/2  in 
magnitude  acting  at  x  =  a/2.   If  a  local  coordinate  system  with  axes 
j,  k,    &  is  assigned  to  the  point  x,  the  force  P/2  will  produce  a 
moment  about  the  &  axis  equal  to: 


m  ;,  = 


(  a/2  -   x)B 


(2-1) 


x=-a/2 


>  X,/ 


P/2  P/2 

Figure  (2-2)  Model  I  showing  arbitrary  point  x 

with  force  P/2  acting  to  the  right  a  distance  a/2  -  x 


Where  B  again  is  the  angle  of  twist  indicated  in  figure  (1-2).   It 


should  be  noted  that  B  is  assumed  to  be  sufficiently  small  so  that 
small  angle  approximations  for  sine  and  cosine  are  applicable. 


The  expression  for 


A2 

d  v 


dx' 


d2v 

dx2 


then  is  (reference  1) 


M, 


EI 
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For  B  the  simple  choice  B  =  Bo  cos 


TTX 


will  be  made.   The  cosine 


expression  is  selected  because  it  is  assumed  that  the  stiffener  will 
be  mounted  such  that  it  is  restrained  from  twisting  at  the  ends  x  = 
a/2  and  x  =  -a/2.      This  is  a  reasonable  assumption  in  the  case  of 
stiffeners  fixed  between  decks  or  frames  of  a  ship. 


d2v 
Utilyzing  the  appropriate  terms  for  M^  and  B,  = —  becomes 

dx 


-.2 

d  y 

dx 


2^y    (  2  ~        COS 
z 


TTX 


(2-2) 


Inserting  this  expression  into  the  formula  for  strain  energy  U  yields 


U 


-a/2    r-p2Bo2 

4EI 


2  ttx   ,  a     ,2  ,  ^^-Bott.2     2  ttx 

cos   (  ■=    -   x)   +  G J  ( )   sin   

a     2  a  a 


,  t,^       t,     2  .  tt  .  4     2  TTX 

+  ECw  Bo   (  -  )   cos   

a  a 


dx 


Carrying  out  the  integration  it  becomes: 

2   2  2 

P  Bo     3,6  +  tt   .  .  __  ,  Bo~   2  a    _^   _  2  .  tt  4  a 

U  =  a   tt-   )  +  GJ  (  — - — )  7  +  ECw  Bo  (  -)  7 

2              a  4  a  4 


EI 


192^ 


2.2   Development  of  the  Virtual  Work  Equation 


Work  is  defined  as  the  movement  of  a  force  through  a  distance. 
In  the  case  of  Model  I  then  the  virtual  work  done  is  obtained  by 
multiplying  the  force  P  times  the  summation  of  all  the  components  of 
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motion  in  the  Z  direction  for  elements  of  the  beam  between  x  =  0  and  x 
=  a/2.   It  is  not  necessary  to  sum  vertical  components  over  the  entire 
length  of  the  beam,  x  =  -a/2  to  x  =  a/2  ,  because  each  half  moves  the 
same  total  distance.   The  expression  for  the  vertical  distance  moved 
by  the  centroid  of  a  cross  section  of  the  beam  at  the  arbitrary  point 
X  is  ( reference  1) : 


A2 

_   d  v     ,  a 

~7~2~      (  2  "  X 
dx 


and  the  work  equation  then  becomes: 

a/2 

W       Vv       B  2_Z_  (  |  _  x  )  dx 

dx 


This  equation  would  suffice  alone  if  the  force  P  were  applied  at  the 
centroid  of  the  cross  section  at  X  =  0;  but  since  we  are  considering 
the  more  realistic  case  where  the  load  is  applied  at  the  outside 
surface  of  the  beam,  a  distance  h/2  above  the  centroid,  the  expression 
must  be  slightly  ammended.   Due  to  the  twist  of  the  section  at    X=0, 
the  load  will  drop  through  an  additional  distance  d  =  h/2(l  -  cos  B^.o  ) 
as  shown  in  figure  (2-3). 
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s~ 

\ 
\ 
\ 
\ 

\B 

\ 

\ 

T 

Figure  (2-3)   Additional  lowering  of  P  due  to  load 
being  applied  a  distance  h/2  above  centroid 

By  the  half  angle  relation, 

2 
1  -  cos  BA=0=  2  sin   B.V2 

~  2(B»V2)2 


Bio/2 


Therefore  d  becomes  approximately 


d  - (h/2)  (Bt0/2) 


(h/4)  Bo2  cos20 


=  (h/4)  Bo' 
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This  term  is  now  multiplied  by  P  and  added  to  the  work  equation  to 
give  a  more  accurate  expression  for  the  work  done. 

W   =    P     O  B-^-2-     (    |— x)dx      +       -^Bo2  (2-3) 


2 

_.  ttx  ,    d    v  w  a        .  ,  Ph    „    2 

Bo    cos    ( 7T  )  (-zr-x)  ax   +    Bo 

dx  4 


2 

d   v                P  ,  a                 .       _                 , 

Again          — 2~    =    2EI —  (~2—    "    x    )       B      '    and 
dx                      z 

when    applied    the    work  equation    becomes: 

?       7  a/2 

P    Bo  C                     2      ttx,     a               .2,           PhD2 

W    =    \  cos         (    -=-    -    x    )       dx    + Bo 

2EIz 
Integrating    this    expression    gives: 


P2Bo2  3        .    6    +    tt2  .  Ph      _    2 

W   =    a  (    —  )       +     — j-    Bo 

EI  96-2 

z 


2.3      Determination    of    the    Critical   Buckling    Load 


To    determine    the    value   of    P    that   will   cause    failure    of    the 
stiffener,    the      critical    buckling    load    Per,    the    principle    of 
conservation   of    energy   will   be    employed.      This    concept    states    that 
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when  a  member  is  gradually  loaded,  the  kinetic  energy  is  zero  and  the 
work  done  by  external  forces  is  equal  to  the  strain  energy  (reference 
2)  : 

6W  =  <5U 


To  determine  6W  and  5U  calculus  of  variations  with  respect  to  Bo  is 
applied  to  the  work  and  strain  energy  equations.  First  to  the  work 
equation : 

2 


W(Bo    +    6 Bo)     =     (Bo    +    6 Bo) 


2    3  2 

Pa       ,    6    +    tt    .     +    Ph_ 


EI 


=    W(Bo    )     +    2BoxBo 


:  96tt  

"23  2 

Pa       ,6_+_L \    +    Ph_ 


EI 


96ti 


5W    =    2Bo6Bo 


rP2a3 

EI 


2       ;  4 

96- 


And    in    like   manner    to    the    strain    energy    equation: 

?       P    ^  fi  +  T  tt     2    a  tt      2    a 

^T      (  2    )+    GJ(a}       I    +    ECw(T}       ! 


z       192 
2    3  2  2  — 

)2+Bo^Bo  |    |y^—     (6    +    2    )+    4^     (GJ    +    ECw(l)2 


Ci 


<5U=2Bo6Bo 


192t 
2 


pV(6±11_)+    £     (gj    +   ECw    (|,2 


LfIz     V192tt2 


4a 


] 


The  two  variational  expressions  may  now  be  equated,  and  the  formula 
for  the  critical  buckling  load  derived. 


6W  =  SU 
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2Bo6Bo 


2  3       2 
Per  a   ,  6  +  ^        Pcrh 


EI       '     rsr       2 


2Bo6Bo 


J""    2  3 

Per  a   ,  6+ 

L  EIz      19 

D 


2tt 


+  EC„(  =  ,2] 


2    3  2 

Pcra  ,6  +  tt       .  Pcrh          tt  "         „„           „       .  tt  .  2   , 

( 2}     +    ~~1 4a"  GJ    +    ECw     (  — ) 

EI                 192tt  4  4a 

z 


■] 


This   quadratic    may    be    solved    for    Per,    and    the    result    expressed    in    the 
form: 


-    (    ^-V) 


6  +  TT 


2    '    a 


iV- 


EI  , 

— S  +    i 

GJ  2 


<  ^4>2  <t)2  EI* 


6  +  TT 


a'  GJ 


+     ( 


192 


(    1    + 


ECw 


6+TT 


GJ 


<  a  r  ) 


1/2      7    /EI    GJ 


(2-4) 


This    formula    corresponds    directly    with    equation    (6-18)    of    reference    1 
for    similar    geometry    and    loading    which,    in    the    notation    of    this 
thesis,    is: 


Per 


y^rcj- 


The  expression  within  braces  in  equation  (2-4)  equates  to  the  non- 
dimensional  factor  IT  of  (2-5).   Figure  (2-4)  shows  the  close  agreement 
between  the  values  of  Y  calculated  from  equation  (2-4)  with  those 
listed  in  tabular  form  in  Table  6-5  of  reference  1. 
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Figure  (2-4)   Agreement  of  calculated  values  of  V  from  eq  (2-4) 
with  tabular  values  listed  in  reference  (1) 
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CHAPTER  3 
MODEL  II  -  TORSIONAL  BUCKLIN3  UNDER  POINT  LATERAL  LOADING 


3.1   Development  of  the  Strain  Enerqy  Equation 


Model  II  again  was  the  flat  bar  tvpe  stiffener  illustrated  in 
figure  (1-4).   In  this  chapter  the  same  type  of  loading  situation  as 
that  of  chapter  2,  a  point  load  coincident  with  the  Z  axis  and  applied 
at  the  upper  surface  of  the  stiffener,  will  be  examined.   A  pictorial 
representation  of  the  loading  will  appear  identical  to  that  shown  in 
figure  (2 -1 ) . 

The  equation  for  the  strain  energy  of  a  flat  bar  in  lateral 
bending  and  torsion  can  be  obtained  in  exactly  the  same  manner  as  that 
for  the  I  beam.   The  only  difference  is  that  the  warping  riqiditv 
term,  Cw,  vanishes.   This  is  because  the  equation  for  Cw  of  an  I  beam 
is: 


Cw  = 


2  3 
tfli  b 
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and  the  flat  bar  stiffener  can  be  likened  to  an  I  beam  with  the  width 
b  of  the  top  and  bottom  flanges  equal  to  zero.   The  strain  enerqy 
expression  then  is: 

a/2, —  ?         ?  2 


U    = 


dx 


dx 


o 


TTX 


The    simple    choice    for    B,    Bo   cos    — —      ,    is    again    made    and    the 


expression    for 


dx2 


is  then  identical  to  (2-1) 


dfv 

dx2 


P     .  a     x  _        TTX 
(  7T  -  x)  Bo  cos  — 
z.  a 


2EI 


The  complete  strain  energy  equation  is: 


a/2 


U   = 


o 


r*^2„    2  2  ,    a  v  2         __    ,    Botk  2       .     2    ttx 

P    Bo         cos      ttx       (    T   -    x)       +    GJ     (    )       sin      —      |     dx 

—  z  a  a 

4EI  a 

i—         z 


G 


which    when    integrated    becomes 


2       2  2  2 

P    Bo            3     ,    6    +    tt       .                  ^T     ,Bott,  a 

U    =      a       ( t—    )       +         GJ     ( )  — 


4EI 


192tt' 
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3.2   Development  of  the  Virtual  Work  Equation 


The  work  equation  is  found  again  by  defining  the  work  done  as  the 
product  of  the  load  P  times  the  summation  of  all  vertical  components 
of  motion  over  one  half  of  the  bar.   The  expression  for  the  vertical 
distance  moved  by  the  centroid  at  any  point  along  the  X  axis  is 
identical  to  the  I  beam  case: 


A2 

d  v    .a     , 

B  —~2       (2  "  X) 
d  x 


Also,  a  correction  term  must  once  more  be  added  to  the  work  equation 

to  account  for  the  load  being  applied  at  the  upper  surface  of  the 

plate  rather  than  at  the  centroid.   The  ultimate  result  is  that  the 

work  equations  for  model  I  and  model  II  are  identical: 

a/2 

P  Bo      \                2  7TX    ,a     ,2  ,     Ph  _  2 
W  =  - — \   cos   (—  -  x)   dx  +  —  Bo 

2EI 
z 

P2Bo2     3   ,6  +JT2.    Ph  _  2 
=  a    ( ~— )  +  —  Bo 

EI  96tt        4 


3.3   Determination  of  the  Critical  Buckling  Load 


By  setting  the  variation  of  the  work  done  equal  to  the  variation 
of  the  strain  energy  the  value  of  Per  is  once  again  determined.   Since 
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the  work  expression  for  the  flat  plate  is  identical  to  that  of  the  I 
beam,  so  too  is  the  equation  for  the  variation  of  the  work: 


2    3  2 

Pa             .6    +    ir    .  Ph 

6W   =    2Bo6Bo  ( y~  )    +    — 

EI  96tt  4 

z 


[ 


1 


While    that    for  6 U    is    slightly    different    owing    to    the    abscence    of    Cw 


U(Bo    +    6  Bo)     =     (Bo    +    6 Bo) 


P!id       {S-LlL)     +    GJ     (1)2    a 
EI  192tt  a         4 


1 


c 


D2    3  2  2 

^-a—       (  %)     +    GJ     (-)       | 

EI  192tt  a         4 

z 


] 


Equating    the    work    and    strain    energy    and    solving    for    Per    yields: 


2Bo6Bo 


6W   =    6U 

2    3  2 

Per    a  /6_+_l_  Ell 


EI 


96tt' 


=    2Bo5Bo 


£££^      (L+-4)    +     gj    (J)'      f 
EI  192tt  a 


PcrV       (6_4_4)     _      £h      +      II!      GJ    =    0 


EI  192tt' 


4a 


When    solved    this   quadratic    gives    for    Per    an    expression    which    may    be 
written    in    a    form    similar    to    equation    (2-2): 
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Per 


24tt' 


6  +  tt' 


o  & 


EI 


GJ 


48tt 


6  +  tt' 


2^2   EI 

■)     (&)      -£ 

a     GJ 


+   ( 


192TT 

6  +  tt' 


1/2 


^  EI  GJ 


(3-1) 


For  typical  stiffeners  the  quantity  ( 


192TT 


)  will  be  several 


2    6  +     FT 

orders  of  magnitude  greater  than  ( — — — T)   (—  r— — — Z. 

6    +    itz  a         GJ 


Therefore,    for 


purposes    of    a    design    formula    it    is    sufficiently    accurate    to    ignore 
this    latter    quantity    and    rewrite    equation    (3-1)     as: 


Per 


I    (192*4       1/2.      _2* 


6    +    it' 


6    + 


a  GJ 


EI    GJ 
z 


(3-2) 


The  formula  for  Per  for  a  flat  plate  stiffener  subjected  to  a 
point  load  in  the  center  of  the  span  is  given  in  reference  (1)  by 
equation  (6-37),  which  in  this  thesis'  notation  is: 


Per  =   16.94 


>/! 


I  GJ 


(1  - 


.87  £ 
a 


EI 


GJ 


Now  to  two-decimal  accuracy  equation  (3-2)  may  be  rewritten  as 


Pc 


r  =   17.16  V- 


EI  GJ 
z 


1  -  .87  £  V_ 


EI 


GJ 


The  formulas  are  in  agreement  within  1.3%. 
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CHAPTER  4 


MODELS  I  AND  II  -  TORSIONAL  BUCKLIN3  UNDER  COMBINED 
LATERAL  AND  AXIAL  POINT  LOADING 


4.1   Description  of  Loading  Condition 


In  chapters  2  and  3  the  failure  by  torsional  buckling  of  a  beam 
subject  to  a  lateral  load  only  has  been  studied;  now  suppose  that  an 
additional  axial  load  is  added  at  the  ends  of  the  beam  as  shown  in 
f  igure  (4-1) . 


Q 


Q 


Figure  (4-1)   Combined  Lateral  and  Axial  Loading  of  a 
Stif f ener 
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The  stiffener  is  again  simply  supported  and  laterally  loaded  by  a 
force  P  applied  at  the  upper  surface  in  the  center  of  the  span,  but  an 
idditional  force  Q  acts  at  the  centroid  of  each  end  to  compress  the 
stiffener  inward.   The  two  forces  are  acting  at  right  angles  to  each 
>ther;  P  along  the  Z  axis  and  Q  along  the  X  axis.   This  dual  axis  type 
oading  would  seem  to  be  much  more  representative  of  the  kind  of 
oading  conditions  to  be  encountered  by  the  frames  of  a  ship.   For 
xample,  two  commonly  encountered  conditions  by  a  ship  at  sea  in  large 
aves  are  hog  and  sag,  shown  in  figure  (4-2). 


HOGGING 


SAGGING 


Figure  (4-2)   Hogging  and  Sagging  Congitions  in 
Large  Waves 
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Here  the  ship  is  supported  by  the  crest  of  a  wave  either  amidships 
(hogging)  or  at  the  ends  (sagging).   In  both  cases  the  vertical 
transverse  frames  of  the  ship,  a  portion  of  which  are  shown  as  dashed 
Lines  in  figure  (4-2),  will  be  placed  in  compression,  i.e.  a  Q  force 
applied.   At  the  same  time  then,  the  sides  of  the  ship  will  very 

ikely  be  receiving  a  lateral  force  P  from  other  waves  in  a  confused 
sea  or  in  a  worst  case  from  collision  with  another  ship,  floating 
object,  or  possibly  from  within  the  ship  from  unsecured  cargo 

uf feting  about.   How  do  the  two  forces  interact,  and  what  effect  does 
the  addition  of  an  axial  force  have  upon  the  capability  of  the 
istiffener  to  accept  a  lateral  load?   These  are  the  questions  this 
:hesis  will  attempt  to  answer  utilizing  the  strain  energy  approach 
/hich  was  used  in  the  previous  chapters  for  a  lateral  load  alone. 


.2   Development  of  the  Strain  Energy  Equation 


In  solving  for  the  strain  energy  of  a  stiffener  subjected  to  both 
ateral  and  axial  loads  the  case  of  model  I,  an  I  beam,  will  be 
:onsidered  first;  so  the  energy  equation  will  contain  a  warping 


"z<^'2  +  -'^>2  +  ECw<— !'2 


dx 


dx 


dx' 
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Again    for    B,    the    angle    of    twist,    the    choice   Bo   cos  H2L    is    made.      The 

a 


expression  for 


dx2 


will  however,  no  longer  be  the  same.   Taking  again 


an  arbitrary  point  x  as  shown  in  figure  (2-2),  there  will  still  be 
present  to  the  right  of  x  the  force  P/2  creating  a  moment  about  the 
local  &  axis  equal  to  the  expression  in  equation  (2-1): 


m   =  T  (  I  "  x)  B 


(2-1) 


There  is  now  though,  another  force  to  the  right  of  x  adding  to  this 
moment  about'-  ,  the  axial  force  Q.   Looking  down  upon  the  stiffener 
from  above  it  is  seen  that  this  addition  to  M  I   is  equal  to  the  force  Q 
times  the  lever  arm  -v  as   shown  in  figure  (4-3): 


Q 


Figure  (4-3)   Additional  moment  Qv  provided  by 
axial  loading 


The  total  moment  M  1   then  is  now 


MS.   =   |  (|  -  x)  B  +  Qv 
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d2 

ind  the  expression  for  — ^—  becomes: 

dx^ 


d  v  MJ 


,  2  EI 

dx  z 


(  ~  -  x  )  Bo  cos  +  —7—         (4-1) 


2EI     v  2      '         a      EI 
z  z 

tfe  now  have  a  linear  second  order  differential  equation  which  must  be 
solved  for  v  in  order  that  the  strain  energy  may  be  found.   We  are 
Looking  for  a  particular  solution  to  this  differential  equation,  and 
:hus  it  will  be  of  the  form  (reference  3): 


Vp  =  C,  cos  — ■ +  c_  sin  — +  C-,  x  cos  ■  ■+  C.  x  sin  (4-2) 

Ia2a3  a4  a 


differentiating    this    expression    for    vptwice   with    respect    to    x    yields: 

£&—      =    -    (2C,    -    -    CJcos    -22      -    H     (2C,    +    2    C9)     sin    -^ 
.,2  a  4        a      1  a  a  3        a      2  a 

dx 

-     (-)       C,    x  cos    — ^-  -     (2)       c.  x     sin      ^X  (4-3) 

a  j  a  a  4  a 


[f    v    is    a    solution    to    equation    (4-1)    the    following    is    true 
P 


,2 

d    Vp              P  ,a              .      _                       TTX     ,        Q            ,_                    TTX        ,      _                    ttx 

7^-  =  -*== —      (=    -    x)    Bo    cos    +  — (c,    cos    +    C„    sin 

-,    2            2EI  2                                      a            EI               1                 a               2                 a 

dx                      z  z 


+    C       x  cos       nX      +    C.     x  sin      VX      )  (4-4) 

3  a  4  a 
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By    equating    coefficients    of    the    trigonometric    functions    in    equations 
(4-3)    and     (4-4),    four   equations    can   be    found    to   be    solved    for    the    four 
constants    C,     through   C..       They   are: 


ci  ■ 


-a   PBo 


4a2Q    +    4EI    t2 
z 


-a   EI    ttpbo 
z 


2  2    2 

(a    Q    +    EI    ttz    Z 
z 


C 


a2PBo 


3  2  2 

2a    Q    +    2EI    tt 
z 


C4    =    ° 


Inserting    these    constants    into   equations     (4-2)     and     (4-4)    gives    for    v 
and    its    second   derivative: 


v      =     ( 


-a    PBo 


4a    Q    +    4EI    tt 
z 


.  TTX 

77    )  cos    — • 
2  a 


,         a    EI    t-PBo 

( z 

2  2    2 

(aZQ    +    EI    tt    ) 
z 


)  sin   — 
a 


+     ( 


a2PBo 


2a    Q    +    2EI    tt 
z 


>  TTX 

tt     )      X     COS     — 

z  a 


(4-5) 


and 


A2 

d   v 


dx' 


P         ,    a  s  _  irx  Q 

=       2EI        (    2       "    X)B°COS-¥     +     EI- 
z  z 


-a  PBo  ,         tt  x 

~~2 2   )cOS  T 

4a  Q  +  4EI    / 

z 


a  EI   nPBo 

-    (  — 5 5-5-)  sin  -— -       +      (- 

(a2Q  +  EIz/)2  a 


2 

a  PBo            ,  tt  x 

— 5— )x  cos  — 

z.  a 


2a  Q  +  2EI  tt 
z 


(4-6) 


d2v 


dx 


PBO  TTX 

cos  


EI 


)  +  (- 


a2Q 


2a  Q  +  2EI   tt 
z 


B 


2 }    "  2   )    X 


3 

a  PBoQ 


(a"Q  +  EIz  tt  ) 


TT  X 

^-rr     sin 

2, 2  a 


(4-7) 


For    simplicity    the    coefficients    in    equation    (4-7)     shall    be    labeled 
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,B,  and  C  as  indicated,  and  (4-7)  then  becomes: 
PBo 


A2 

d  v 


dx 


EI 


cos  (A  +  Bx)  -  C  sin  


nd 


(  dx2   ) 


(  p°_  }  2  cos2  _2L  (  A2  +  2ABx  +  B2  x  2 

Li  a 

z 


0   PBo   ,.^  ,  _,_  »       ttx     2  tx 
-  2  -rr: (AC  +  BCx)  sin  cos  


EI 


„2   .  2  -x 

+   C   sin   

a 


(4-8) 


ow  with  equation  (4-8)  there  is  sufficient  information  to  solve  the 


train  energy  equation. 
ra/2  ' 
U  -\ 


J  o 

■a/2 


2  -   2 

r^-r   /  d  v  N2    __,dB  <2    __   .   d  B.  2 

EI   (  s  )   +  GJ(^ )  +  ECw  (  ~-  ) 

z    ,  2.  dx  j.. 2 


dx 


dx' 


dx 


r~   2n  2 
p  Bo 


o   L_ 


EI 


cos2  — —   (  A2  +  2ABx  +  B2x2) 


-  2PBo  (  AC  +  BCx  )  sin 


■X  TTX 

—  cos  


^T  ^2   .2  ~x     __  ,  Bo1"^      ti  : 

+  EI  C   sin  +  GJ  ( )   sin  

z         a  a  a 


+  ECwBo   (  -  )   cos   

a  a 


dx 


(4-9) 


arrying  out  the  integration  of  equation  (4-9)  and  reinserting  the 
roper  quantities  for  A,B,  and  C  gives  as  a  final  expression  for 
train  energy: 
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u  = 


2   2  3 
P  Bo  a 

EI 


(6  +  \    )     -     (6  +  \     )  ( 


a2Q 


192 


24tt2    4a2Q  +  4EI  r2 

z 


)  + 


2        2 

(i_±JL_)  (_   a  Q 


12u2     4a2Q  +  4EI  it2 

z 


2   2  5 
P  Bo  a 


(- 


2Q 


2  2  2 

(4a  Q  +  4EI  tt  ) 
z 


2   2  7 
P  Bo  a 


2  2 
64EI  tt    Q 
z 


2  2  4 

4a  Q  +  4EI  t  ) 4 
z 


8a2Q 


2  2  3 

(4a  Q  +  4EI  O  J 
z 


+     GJ(- 


BoT 


+    ECwBo2 (-) 4 (f ) 
a    4 


(4-10) 


4.3   Development  of  the  Work  Equation 


The  work  done  will  again  be  found  from  the  principle  that  it  is 

force  times  distance.   The  work  performed  by  the  load  P  is  given  by 

the  same  expression  used  in  chapter  2: 

a/2 


W  =  P 


D- 


A2 

d  v 


dx 


(j   -  x)dx    +   -^-  Bo 


(2-3) 


Now,  however,  there  will  be  a  new  expression  for 


A2 

d  v 

A       2 

dx 


,  that  of 


equation  (4-6) .   Also,  there  must  now  also  be  added  the  work  of  the 
load  Q  at  the  end  x  =  a/2  as  it  moves  through  a  distance  -u   (figure 
4-4)  . 
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Figure  (4-4)   Motion  of  force  Q  through  distance  -u 


The  total  work  done  by  both  loads  will  be: 
a/2 


W  =  P 


o  d  v   ,a     .  ,      Ph   „  2    _ 

B  —j-  {-2    -  x)dx  +  -^  Bo   +  Qux=a/2 

dx 


(4-11) 


All  the  necessary  quantitites  in  equation  (4-11)  are  known  with  the 

exception  of  u    ,  -   To  find  this  unknown  the  definition  of  membrane 
x=a/2 

strain  (reference  5)  will  be  used: 

-du 


dx 


+  I    (  _^L  )2  +  I    (-dw  )    2 
2    {      dx    '  2  v  dx  ; 


Utilyzing  inextensibility  assumptions  this  membrane  strain  is  set  to 
zero: 


_  _  n     -du     1  ,  dv  .2    1  .    dw   .  2 
C  "  °  =  "dx—  +  2  ("dx-}   +  2  (^x~} 


or 


du    _  1  .    dv   .  2    1  ,   dw   .  2 
dx    "  2  K    dx   j      2  ^   dx   ; 


(4-12) 
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As  we  are  dealing  with  failure  of  the  stiffener  by  a  torsional 
(twisting)  mode  augmented  by  an  axial  compressive  force,  it  is 
expected  that  deflection  in  the  Z  direction  will  be  minor  compared 
with  that  in  the  Y  direction.   Therefore,  for  purposes  of  a  design 
relationship  it  will  be  considered  sufficiently  accurate  to  ignore  the 
— j term  in  equation  (4-12) ,  and  the  relation  becomes: 

du      1     dv_   2  3) 


dx      2      dx 

An  expression  for  u  may  now  be  found  by  solving  equation  (4-13).   Both 
sides  are  multiplied  by  dx  and  then  integrated: 


/i,,   -    1  /   dv   n  2   , 

du   -     2  (  "dir }     dx 


To  find  the  total  work  done  by  Q,  it  is  necessary  to  find  the 
deflection  of  the  extreme  right  hand  side  of  the  beam,  at  x=a/2.   So 
the  limits  of  integration  must  be  over  the  entire  length  of  the  beam, 
x=-a/2  to  x=a/2;  or,  twice  over  half  the  beam,  x=0  to  x=a/2: 

a/2 


For  v  equation  (4-5)  is  used,  and  the  first  derivative  of  v  with 
respect  to  x  is: 
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dv 


dx 


=    ( 


2™ 

Tra    PBo 


4a2Q   +    4EI    tt2 
z 


)    sin 


2  2 

a  EI    tt    PBo 

—     "    (      9      +    2 *-*>    cos    -El 

a  (a2Q   +   EVV 


D 


E 


)      COS 


TTX 


aTrpBo 


)    x   sin 


TTX 


(4-35) 


\gain    for    simplicity    the    coefficients    in    equation    (4-15)    are    labeled 
D,E,F,    and  G    as    shown;    and    equation    (4-14)    becomes: 


.a/2        r 


u 


x=a/2 


.     2         TTx       /r^2         __,  Jl    2, 

sin        (D      -    2DGx    +    G    x    ) 


+    2D(F-E)     sin    — —cos    -^-    +    2G(E-F)    x    sin    — —    cos  --^ 


+    cos2    — ^—      (    F2    -    2EF    +    E2) 


dx 


(4-16) 


Carrying    out    the    integration    of    equation    (4-16)     and    replacing 


coefficients   gives    for    the    displacement    u    _        : 

x — a/  c. 


5„2„    2 

a    P    Bo 


u^_ 


x=a/2 


(4a2Q    +    4EI    v2) 
z 


1         12         ; 


24EI    tt' 
z 


2  2 

4a    Q    +    4EI    tt 
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+     ( 


8EI    tt' 
z 


2  2 

4a    Q    +    4EI       tt 
z 


(4-17) 
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All  quantities  are  now  known  to  determine  the  work  done  by  the 


combined  lateral  and  axial  loads.   Utilyzing  the  proper  expressions 

d2 
for  B,  — ^   (eq.  4-6),  and  u       (eq.  4-17),  equation  (4-11)  becomes: 

dx2  x=a/2 


a/2 


W=  P 


ttx   3    P       ,  a      v  _       TT  x 
Bo  cos  — : —  <  -tt== —    (  ~  -  x  )  Bo  cos  — — 
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EI 
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4a  Q  +  4EI  tt 
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TTX 
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Integrating  equation  (4-18)  yields  for  work: 
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4.4   Determination  of  the  Critical  Buckling  Load  Equation 


To  determine  the  relationship  between  the  loads  P  and  Q  and  how 
they  will  affect  buckling  of  a  stiffener,  the  conservation  of  energy 
principle  will  once  more  be  employed.   First  calculus  of  variations 
with  respect  to  Bo  is  applied  to  equations  (4-19)  and  (4-10)  to 
determine  as  before  6W  and  5U: 
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By  conservation  of  energy  6 W  =  6U;  thus  equations  (4-20)  and  (4-21) 
may  be  equated,  and  the  ultimate  relation  between  P  and  Q  appears  as 
the  following  quadratic  with  respect  to  P: 


o2  3    c     2 

P  a  .    6  +  7T   . 
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Equation  (4-22)  presents  the  relationship  between  P  and  Q  when  they 

are  applied  to  an  I  beam.   If  both  forces  are  applied  in  like  manner 

to  model  II,  the  flat  bar  stiffener,  the  only  difference  in  equation 

4 

TT 

(4-22)  is  that  the  final  term,  ECw ( -) ,  will  be  absent  since  Cw 


for  the  flat  bar  is  zero. 
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CHAPTER  5 
EVALUATION  OF  RESULTS  OF  THE 
BUCKLING  LOAD  EQUATION 


5.1   Nond imensionalization  of  the  Buckling  Load  Equation 


It  is  now  our  desire  to  use  equation  (4-22)  to  study  the 
interaction  between  the  lateral  and  axial  loads  P  and  Q  and  determine 
how  changes  in  the  physical  dimensions  of  a  stiffener  effect  its 
ability  to  withstand  against  failure.   The  most  interesting  comparison 
to  be  made  from  equation  (4-22)  would  be  a  plot  of  lateral  load  P 
versus  axial  load  Q,  i.e.  a  means  to  predict  for  a  given  amount  of 
lateral  load  how  much  simultaneous  axial  load  a  stiffener  will  take 
before  buckling  and  vice  versa. 

In  order  to  construct  such  a  plot  it  will  be  convenient  to  first 
nondimensicnalize  equation  (4-22).   An  ideal  means  to 
nondimensionalize  P  and  Q  is  to  express  them  as  a  fraction  of  the 
Euler  buckling  load  (equation  2-5  of  reference  1) : 

TT   2EI 

n   -        z 

Qe 


2 
a 
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Qe 


The  Euler  critical  buckling  load  Qe  is  that  load  applied  as  in  the 

manner  shown  in  figure  5-1  that  will  just  cause  in-plane  deflection  of 

the  beam  without  the  twisting  effects  of  torsion. 

.-Z 


->  X 


Figure  (5-1)  Euler  Buckling  of  a  Beam 


Use  of  the  Euler  buckling  load  as  a  nondimensionalizing  parameter  is 
particularly  well  suited  to  equation  (4-22)  noticing  the  frequent 
occurrence  in  it  of  the  expression 

2  2 

4a  Q  +  4EI  tt 


which  may  be  written,  if  Q  is  taken  as  positive  in  tension  as: 

Q   » 


2       2        2 
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The    final    dimensionless    form   of    (4-22)    will    be: 
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(5-1) 


where 


L  =    1  - 


Qe 


To  three  decimal  accuracy  equation  (5-1)  becomes: 


(-V 

Qe; 


16    - 


Q 


Qe 


1.433     (h2      +      3.701     (h3    -    2.467(i)4 

i-i  J_i  Li 


♦  <&2 


13  14 

1.234     (±)J    +    2.467     (±) 


+    2.467  ( 


P    ,     h 


949    ~    "    24.352     ( — ~ — )    =      0 
z  a    I 


Qe         a 


(5-2) 


Equation  (5-2)  will  be  applicable  to  an  I  beam,  and  to  make  it  useful 
for  a  flat  bar  stiffener  as  usual  set  Cw  to  zero  and  delete  the  final 
term. 


5.2   Model  I  -  Effects  of  Varying  Stiffener  Dimensions  on 
the  Value  of  the  Critical  Buckling  Load 


The  results  of  equation  (5-2)  for  I  beams  of  various  sizes  are 
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shown  graphically  in  figures  (5-2)  to  (5-5) .   The  points  on  the  graphs 
were  derived  from  computer  analysis  shown  in  Appendix  A.   Along  the 
abscissa  of  each  of  the  graphs  is  shown  the  dimensionless  ratio  of 
axial   load  Q  to  the  Euler  buckling  load  Qe,  while  the  ordinate  gives 
the  ratio  of  the  critical  torsional  buckling  load  Per  to  the  Euler 
load.   Thus  by  calculating  the  Euler  buckling  load  of  a  stiffener  and 
being  given  a  value  of  axial  load  Q,  we  can  enter  the  graphs  with  a 
value  of  Q/Qe  and  obtain  a  prediction  of  how  much  lateral  force  Per, 
applied  at  the  center  of  the  span,  that  the  stiffener  will  withstand 
before  torsional  buckling  occurs.   For  example,  in  figure  (5-2)  if  an 
axial  force  equal  to  a  value  of  ,4Qe  is  applied  to  stiffener  3  with 
h/a  ratio  equal  to  .14,  then  a  lateral  force  of  just  .16  Qe  should 
produce  failure  of  the  stiffener  by  torsional  buckling. 

Though  no  previous  design  formulas  relating  simultaneous  axial 
and  lateral  loading  of  a  stiffener  to  cause  torsional  buckling  were 
found,  the  results  of  equation  (5-2)  do  appear  reasonable.   It  is 
possible  to  verify  the  predicted  results  of  Pcr/Qe  at  the  intersection 
with  the  ordinate  (Q/Qe=0).   For  example,  again  for  stiffener  3  of 
figure  (5-2)  with  h/a  equal  to  .14,  the  complete  set  of  scantlings  for 
this  I  beam  are: 


a=3.07m 

h=.43m 

b=.09m 

tf=.02m 

tw=.007m 
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1.0- 


Pcr/Qe 


Numbers  adjacent  to  each 
curve  correspond  to  stiffener 
numbers  in  Appendix  A 


h/a  increasing 


Figure  (5-2)   Pcr/Qe  vs  O/Qe  for  I  beam  varying  h/a 
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Pcr/Qe 


.0 

.8 

.6 

4 
s/   5 

6 

_.    h/tw   increasing 

.4 

.2 

i 

i               i              r~~ — - — - — ; 

4        .6        .8        1.0 


O/Qe 


Figure  (5-3)   Pcr/Oe  vs  Q/Oe  for  I  beam  varying  h/tw 
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1.0- 


Pcr/Oe 


b/tf  increasing 


1.0 


Q/Qe 


Figure  (5-4)   Pcr/Oe  vs  Q/Qe  for  I  beam  varying  b/tf 
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1.01- 


.8  - 


Pcr/Qe 


tf/tw  increasing 


1.0 


Q/Qe 


Figure  (5-5)   Pcr/Qe  vs  Q/Qe  for  I  beam  varying  tf/tw 


Utilyzing  these  dimensions  the  warping  constant,  St.  Venant's  torsion 
constant,  and  moment  of  inertia  of  the  stiffener  about  the  web  are: 

r„  .    (.02m) (.43m)2(.09m)3     .  . _  v  -.-7   6 
Cw  =  _  .  =I.lzJxlO    m 

2(. 09m) (.02m)3  +  (. 43m) (. 007m) 3    _  OQ„    . n-7   4 

t       2(.09m)3(.02m)     0  ,-,n    , n-6  4 

I    =  t-z L    =      2.430  x  10   m 

z  12 


Now  the  Euler  buckling  load  Qe  for  this  beam  calculated  from  equation 
(2-5)  of  reference  (1)  is: 

7T   EI 

Qe  =  ■ 


2 
a 

tt2(2.068  x  108  KN/m2)   (2.430  x  10~6m4) 

(3.07m)2 

=  526.2  KN 


To  find  the  critical  load  for  torsional  buckling  from  a  lateral  load 
alone  we  use  equation  (6-18)  of  reference  (1): 


ytfEIzGJ 

Per  =   


2 
a 


0         is  obtained  by  entering  figure  (2-4)  of  this  thesis  with 

2    GJ   _   (3.07m) 2 (7.954  x  1Q7  KN/m2)  (5.292  x  10~7m4) 
ECW       (2.068  x  108  KN/m2)  (1.123  x  10_7m6) 

=  17.08 
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and  reading  off  0    of  15.3.   Now  the  critical  torsional  buckling  load 

is: 

pcr  _  15.3  C(2.068xlQ8KN/m2) (2.430xl0"6m4) (7.954xl07KN/m2) (5.292xl0"7m4)  ]  1/2 

(3.07m)2 

=  23  6.1  KN 

A  calculated  ratio  of  Pcr/Qe  then  is: 

Pcr    _    23  6.1  KN 
Qe         526.2  KN 

=  .449 

The  value  of  Pcr/Qe  obtained  from  equation  (5-2)  for  this  same 
stiffener  and  shown  on  figure  (5-2)  for  Q/Qe=0  is  .455;  an  agreement 
of  1.3%.   Table  5-1  lists  the  values  for  Pcr/Qe  at  Q/Qe=0  for  all  the 
curves  of  figures  (5-2)  through  (5-5)  and  shows  the  good  correlation 
between  the  results  of  equation  (5-2)  and  the  values  calculated  using 
equations  (2-5)  and  (6-18)  of  reference  (1).   All  cases  are  within  2% 
agreement  except  stiffener  1 ,  figure  (5-2),  and  stiffener  9,  figure 
(5-4).   In  these  two  cases  a  problem  arises  in  trying  to  pick  off  an 

accurate  value  of  If   from  figure  (2-4)  for  use  in  equation  (6-18).   We 

2 
are  in  the  vicinity  of  a  GJ  equal  to  one  on  figure  (2-4) .   Here  the 

ECw 
values  of  ¥   are  changing  quite  rapidly,  and  a  small  error  in  If  can 

lead  to  an  error  of  several  percentage  points  in  Pcr/Qe. 
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Curve 


Pcr/Qe  at  Q/Qe=0 
equation  (5-2) 


Pcr/Qe  AT  Q/Qe=0 
reference  (1) 


Percent 
Ag  reement 


Fig 

5-2  ,Beam 

1 

.986 

Beam 

2 

.586 

Beam 

3 

.455 

Fig 

5-3 ,Beam 

4 

.596 

Beam 

5 

.586 

Beam 

6 

.573 

Fig 

5-4 ,Beam 

7 

.945 

Beam 

8 

.586 

Beam 

9 

.486 

Fig 

5-5  ,Beam 

10 

.748 

Beam 

11 

.586 

Beam 

12 

.575 

1.056 
.583 
.449 
.598 
.583 
.575 
.938 
.583 
.515 
.734 
.583 
.587 


6.6% 
.5% 

1.3% 
.3% 
.5% 
.5% 
.7% 
.5% 

5.6% 

1.9% 
.5% 

2.0% 


Table    (5-1)       Comparison   of    values    for    Pcr/Qe    at    Q/Qe=0    for    an    I 
beam    as    calculated    by    equation    (5-2)    versus    values    calculated    by 
equations    (2-5)     and     (6-18)    of    reference    (1). 


At    the    right    hand    side    of    the    curves   of    figures    (5-2)    through 
(5-5)    the   graphs   asymptotically    approach    the    x    axis,    but    have    no   value 
at    Q/Qe   =    1.      This    is    because    at    Q/Qe=l    the    expression    1/L    in    equation 
(5-2)    will    have    the    indeterminate    form   1/0.      This    is    to    be    expected 
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though,  because  at  Q/Qe=l  the  stiffener  has  failed  due  to  Euler 
buckling,  and  equation  (5-2)  can  no  longer  be  used  to  describe  the 
relation  between  Per  and  Q. 

The  general  trend  of  all  the  curves  is  similar,  a  rapidly 
decreasing  initial  slope  with  the  values  for  Pcr/Qe  falling  off 
approximately  50%  by  Q/Qe   =.37,  then  tapering  off  of  the  slope  as  the 
curves  approach  Q/Qe=l.   The  d imensionless  parameters  that  were  varied 
in  constructing  the  graphs  were  h/a,  h/tw,  b/tf,  and  tf/tw.   In  each 
figure  the  varied  parameter  has  been  doubled  twice  to  construct  the 
three  curves.   It  is  seen  that  the  ratios  which-  exhibit  the  greatest 
influence  on  Pcr/Qe  are  h/a  and  b/tf.   Varying  h/tw  has  practicallv  no 
effect  upon  Pcr/Qe,  while  changing  tf/tw  at  first  appears  to  have  a 
fairly  significant  effect  but  seems  to  steady  out  after  tf/tw  gets 
above  about  2.8. 


5.3   Model  II  -  Effects  of  Varying  Stiffener  Dimensions  on  the 
Value  of  the  Critical  Buckling  Load 


Figure  (5-6)  shows  the  results  of  equation  (5-2)  for  model  II, 
the  flat  bar  stiffener.   Once  again  the  values  of  Pcr/Qe  for  Q/Qe=0 
may  be  checked  against  those  obtained  using  the  formulas  of 
Timoshenko.   Consider  for  example  stiffener  13  with  h/a=.0523.   The 
stiffener  will  have  dimensions  and  factors  as  follows: 
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Pcr/Qe 


h/a  increasing 


1.0 


Q/Qe 


Figure  (5-6)   Pcr/Qe  vs  Q/Qe  for  Flat  Bar  varying  h/a 
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a  =  3.059m 


h  =  .16m 


tw  =  .0064m 


Cw  =  0 
3   _       (.16.) (.0064m)3     =   1-398  x  1Q-8m4 

z      =  (.16m) (0064m)3   _  -9  m4 

z  12 


The  Euler  buckling  load  from  equation  (2-5)  of  reference  (1)  is 

tt   (2.068  x  108KN/m2)  (3.495  x  10~9m4) 
(3.059m)2 

=  .762KN 


Per  for  the  flat  bar  from  equation  (6-37)  of  reference  (1)  is 

16.94  J  EI  GJ  .    |  EI 

Per  -  J— £-    (1  -  .87  £  v  -4-  ) 

2  a  V   GJ 


=1.56  4  KN 
Per       1.564  KN 


Qe         .762  KN 


=  2.052 


The  value  of  Pcr/Qe  obtained  from  equation  (5-2)  is  2.079,  an 
agreement  of  1.3%.   The  agreement  of  the  values  for  the  other  two 
curves  are  shown  in  table  (5-2)  .   It  is  seen  that  there  is  good 
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agreement  in  all  three  cases. 


Curve  Pcr/Qe  AT  Q/Qe=0     Pcr/Qe  AT  Q/Qe=0     Percent 

equation  (5-2)        reference  (1)     Agreement 

FIG  5-6, Beam  13  2.079  2.051  1.3% 


FIG  5-6, Beam  14  2.005  1.978  1.4% 


FIG  5-6, Beam  15  1.864  1.817  2.6% 


TABLE  (5-2)  Comparison  of  values  of  Pcr/Qe  at  Q/Qe=0  for  a  flat  bar 
as  calculated  by  equation  (5-2)  versus  values  calculated  by  equation 
(6-37)  of  reference  (1). 


As  before  with  the  I  beam,  equation  (5-2)  is  undefined  for  the 
flat  bar  at  Q/Qe=l  owing  to  the  quantity  1/L  being  equal  to  1/0;  but 
this  is  acceptable  since  the  bar  has  undergone  Euler  buckling  at 
Q/Qe=l.   The  shape  and  trends  of  the  curves  are  identical  to  the  I 
beam,  and  again  the  curve  falls  50%  by  about  Q/Qe=.36.   It  is 
interesting  to  note  that  whereas  for  the  I  beam  Pcr/Qe  is  always  less 
than  one,  Pcr/Qe  for  the  flat  plate  starts  out  greater  than  one  by  a 
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factor  of  approximately  two  and  remains  greater  than  one  until 

Q/Qe=. 29.   This  would  indicate  the  much  lower  relative  Qe  for  the  flat 

bar  and  its  much  greater  susceptibility  to  failure  by  Euler  buckling 
than  the  I  beam. 


Owing  to  the  geometric  simplicity  of  the  flat  bar,  the  only 
dimensionless  parameter  that  may  be  varied  independently  is  h/a.   If 
h/tw  is  varied  while  keeping  h/a  constant,  there  will  be  no  change  in 
equation  (5-2) .   As  Cw  is  zero  and  the  ratio  J/Iz  for  a  flat  bar  is 
always  4,  the  only  term  that  will  alter  equation  (5-2)  is  h/a. 
Moreover,  as  can  be  seen  from  figure  (5-6),  even  the  changing  of  h/a 
has  very  little  effect  upon  Pcr/Qe. 
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CHAPTER  6 


CONCLUSIONS 


In  this  thesis  has  been  studied  the  problem  of  torsional  buckling 
of  I  beam  and  flat  bar  type  stiffeners  under  combined  lateral  and 
axial  loads.   First  examined  was  the  case  of  lateral  loading  alone. 
An  energy  method  was  employed  to  determine  critical  buckling  stress 
Per  of  a  stiffener,  and  equations  (2-4)  and  (3-1)  were  developed  to 
cover  lateral  loading  of  I  beams  and  flat  bar  stiffeners  respectively. 
Good  agreement  was  shown  between  equations  (2-4)  and  (3-1)  and  the 
published  formulas  of  Timoshenko.   One  distinct  advantage  of  equation 
(2-4)  for  an  I  beam  is  that  it  gives  an  exact  formula  for  the 
coefficient  X  in  reference  (1)  equation  (6-18)  rather  than  a  tabular 
listing  of  ¥   for  different  sizes  of  stiffeners.   As  was  noted  in 
chapter  5,  when  trying  to  verify  some  of  the  values  of  Pcr/Qe  it 
becomes  quite  difficult  to  pick  off  an  accurate  value  of  tf  when  near 
the  bottom  of  the  curve  in  figure  (2-4)  .   In  particular,  for  the  two 
cases  in  Table  (5-1)  where  the  percent  agreement  between  Pcr/Qe  as 
calculated  from  this  thesis*  equations  and  those  of  reference  (1)  were 
6.6%  and  5.6%,  the  differences  fall  to  just  .5%  and  .2%  when  "X  is 
calculated  from  equation  (2-4)  rather  than  read  from  the  graph. 


62 


Next  was  approached  the  problem  of  combined  lateral  and  axial 
loading  of  a  stiffener.   Again  an  energy  approach  to  the  problem  was 
used,  and  the  ultimate  result  was  equation  (4-22)  which  when 
nondimensionalized  using  Euler  buckling  load  took  the  form  of  equation 
(5-1).   Though  no  other  previous  formula  relating  torsional  buckling 
under  combined  loads  was  found  for  comparison,  the  end  points  of  the 
curves  of  figures  (5-2)  through  (5-6)  could  be  checked  and  were  found 
to  be  in  good  agreement  with  results  obtained  using  other  published 
formulas.   The  trends  of  all  the  curves  appear  reasonable:   a  constant 
decrease,  more  rapid  at  first,  of  the  amount  of  lateral  load  a 
stiffener  can  sustain  as  an  increasing  axial  load  is  applied.   As 
another  vote  of  confidence,  if  Q  is  set  to  zero  in  equation  (4-22)  the 
expression  reverts  to  that  of  equation  (2-4)  which  was  verified 
against  equation  (6-18)  of  reference  (1) .   An  area  for  research  would 
be  to  experimentally  verify  the  results  of  equation  (4-22),  though  it 
would  certainly  be  expensive  due  to  the  complexity  of  a  test  apparatus 
capable  of  imposing  and  measuring  simultaneous  lateral  and  axial 
loads. 


Equation  (4-22)  is  somewhat  lengthy  for  hand  calculations, 
particularly  when  the  constants  Cw,  H,  and  Iz  must  also  be  found 
before  using  it.   It  is  not  too  cumbersome  if  only  one  stiffener  is 
involved,  but  if  a  number  of  stiffeners  are  being  investigated  it  does 
lend  itself  well  to  a  computerized  study,  especially  if  an 
investigation  is  to  be  done  varying  stiffener  dimensions.   The  simple 
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BASIC    program   with    accompanying    results    shown    in  Appendix  A  was 
written    and    utilyzed    to  generate    the   values    for    constructing    figures 
5-2    through    5-6. 
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APPENDIX  A 


Computer  Adaptation  of  Equation  (4-22) 


0  IN!".:  r  N,A,B  ,H,TF  ,   ■ 

1  REM  N  -••'=  ST  IF-!  'ENI     '  IMBER 

2  REM  A  -    STIFF!         :r-  !  ' 

3  REM  "'    =  STIFFENER  FLANGE  BREADTH 

4  pfm  w  =  STIFFENER  HEIGHT 

5  REM  TF  =    STIFF,      FLANGE  THICKNESS 
'  6    R  E  M  ;  W  =  S  "i  I  F  F  E  N  E  R  W  E  B  T !  -i  I  C  K  N  E  S  S 

7  REM  •■  -~0  FOR  E  BE  if  ,  1  (AR 

JO  LET  CW=  (Ti  '*  (H      ■  E<"""'3)  )  ■ 
:0  LET  J=<  (2*E*  (TF'"3)  >  +  (H      '3;  )  )  '   ' 

0  LET  I  =  <  (  < B--3 )  *TF )  /  6 )  +  (  <  <  H*  *  rW""3  >  )  / 1 2  )  *X  ) 

0  FOR  Q=Q  TO  . S  STEP  . 2 

4  LET  L-l-Q 

[0  LET  Ai  -Q*  (  ( 1 ,  433*  (  <  i/L)  "■•2)  )  +  i/L)'v,3))~<    '         ;   "4)  )  5 

10  LET  A2-=  (  ;":  "  2)  *  <  ( 1 ,  234*  (  (  1  /L)  A3)  )  +  (2-  1-67*  •   (J.  /L)  ""h ) 

'0  LET  f'  5=.  B  i  &--A1+A2 

&0  LET  31=        (H 

1 0  L  ET  C=  <      ;=*  ( J  /  I )  )  --  ( 24 .,  352*Cl  I   ( A  ■'■"2)  *  I)  ) 
.20  !  .ET  P  =  «  B 1  ■*•<('<  B 1  "••2 )  -  <  4  *A3  *C )  )  "- .  5 )  )  /  ( 2*  A3 ) 

!§o    !  :  ' 

;|5  ..PRINT  "S  I  CFFENER 
q   LPRIN1  ,:';  IFI  .  !! 

.50   ■  ;  •   i  Ti 

Bf         •      r.'P  "  '  ■  n   : ■'"•-.       np~"  ■   7~.       i!U II   .,   !...!       IIXCTj-n   ..   p        "  Tl'l  "f"l  =  i 

/-      LPT  lN"! 
,|o    LPF  INT     "H    i         " ;  (H/A  '  ;  .CH/TW)  ,  "B/TF     ='   ,  '  '        TF 

!00    LPR  [NT    "   '    '" -'      '.  .   "  ■"    .:' 

IK       1EXT    Q 
BO    END 
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riFT-'E!  i! 
riFFE  JE 
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(DIMENSIONS    AF 
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W=  .007 
4 , 50000 1 


/QE 
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.  4 
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/QE  -~  ) 
•  .,  2 

'  '  ==  .  4 


P/PJ 
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1  5 1 
379R3 
3„  30290* 

:   -/-ri  ?AQP- 
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flFFENER  DIMENSIi 


/A      =     ,  1 400652 

TF/TW    "    2.35714 


}>:•    9o  OOOOC 

N/TW 


TW=     .  < '."■' 


■:•      iOOI 


fUh. 

/QE: 


/QE 


p  /per    =  „  4*34944  1 

P/QE    -  .3430542 

JE    =  .1632 

P/QE    =  5„  78763 IE --02 

P/QE     =  1.         52       ; 
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TIFFJ  !  ILJi  ii   El 


riFR 


NEF     DIMEN  iREs 


B~    9, 


■-).- 


/A    -    ,  231'    "' 
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0 
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H/  TW 


a  .i  „ 


/QE 
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n    -.J  /  JOwf 

p 

/  f  :;i  p 

= 

- 

: 

= 

v7t^9C;=; 

p 

/  f .  s  ■*  ■• 

!! 

. 

= 

1  ,  7082  :  .:.'    ■ 

w 


B/TF 


■ 


.'.;'  XmX)  J 


69 


TIFF  Eh;    ■•     NUMBER     13 
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